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A CLASS OF de LAVAL NOZZLES* 

By S. V. Falkovich 

A study is made herein of the irrotational adiabatic motion of 
a gas in the transition from subsonic to supersonic velocities. A 
shape of the de Laval nozzle is given, which transforms a homoge- 
neous plane-parallel flow at large subsonic velocity into a super- 
sonic flow without any shock waves beyond the transition line from 
the subsonic to the supersonic regions of flow. The method of 
solution is based on integration near the transition line of the 
gas equations of motion in the form investigated by 
S. A. Christianovich (reference 1). 


1. Fundamental equations. - A plane, steady, irrotational, 
adiabatic flow of a gas is considered. In thiB case, the equations 
of motion, as is known, have the form 


dpu , dpv _ A 

sr + aT 


By 3u 

SE " ^ 


= o 


(i.i) 


W 2 A K B‘ K P 0 
2 + K-l p K-l P q 


( 1 . 2 ) 


where 

p density of gas 

u,v components of velocity along x- and y-axes 
p pressure 

W absolute value of velocity 
K adiabatic exponent 


Subscript 0 denotes condition of gas at rest. 


. *"Institut Mekhaniki Akademii Nauk Siuza, SSR" Prikladnaia 
Matematika i Mekhanika, Tom XI, 1947. 
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From equations (l.l) there exist two functions, the velocity 
potential cp (x,y) and the stream function \Kx,y), determined 
by the equations _ 


d<P * u dx + v dy dty = — (-v dx + u dy) (1.3) 

p 0 

If u = ¥ cos 9 and v = ¥ sin 9, where 9 is the angle 
between the velocity vector and the x-axis, are substituted in 
equations (1.3) and are solved for dx and dy, 




w 


P ¥ 


a , . „ + 22 SSS.2 


w 


P w 


(1.4) 


The concepts x and y and <p and are functions of the 
variables ¥ and 9; then 

dqp * d¥ + || d6 


dt = d¥ + d© 

o¥ 09 


If these expressions are substituted for d<P and d'l' in equa- 
tions (1.4), 


dx 


dy 


_ ( ooa 0 dV _ fo sin 9 3\jA ( cos 9 3qp _ fo sin 0 3}jf\ 

\¥^¥ p ¥ 3¥/ \ ¥ 3© p ¥ 3©/ 

„ ( B±a 9 £>? + P_ COB 9 dw + / sin 9 3<P + f_0 cos 9 3^ 

\ w 3¥ Pq ¥ 3 ¥/ \ w 3© p w 3©/ 


(1.5) 

In order for dx and dy, determined frcan equations (1.5), to be 
total differentials, it is necessary end sufficient that the fol- 
lowing equalities hold: 
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a 

( cos 

_e 

a«P _ 

fo 

sin 6 

dt\ 

_ a 


V* 


w 

P 

. V 



a 

( sin 

_e 

&P . 

fo 

cos e 

dt\ 

a 


Vw 


+ 

p 

¥ 

3wy 

“ ^ 

i equation 

(1.2) and the 

condition 




( 1 . 6 ) 


p W E 


1/ 


4. (1 0 
d¥ Ip 


_2 X 

p a 2 


where 

a velocity of sound 

If the differentiation in equations (1.6) is carried out. 


dcp 






p-» 

where M = w/a is the Mach number. 


(1.7) 


In equations (1.7), as the independent variable, in place of 
the velocity ¥ a new variable s (reference l) is introduced, 
which is connected with ¥ by the relation 


s = 


Va 2 -W 2 


dW 


/¥ a ¥ 

Equations (1.7) then assume the form 


a<? aKF** 

~3e 3S 


.rr at 


where 


*0) 


k(b ) “ ( p ) (1 “ M > 


( 1 . 8 ) 


(1.9) 


( 1 . 10 ) 


inasmuch as equation (1.8) is a function of s. 
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The system of equations (1.9) in the regions vhere the flow 
velocity is subsonic, W < a, is of the elliptic type* Hence, 
any solution flp(0,s) and * ^(e,s) of equations (1.9) rep 
resents two functions analytio in the variables s and 0 up to 
the line of transition to the region of supersonic velocities. 

If the Jacobian of the transformation of the region in the 
plane 0,s on the plane q>, ^ 



does not become zero over a certain segment of the transition line 
and therefore also in a certain region on the supersonic side, the 
functions q>(0,s) and >|f(0,s) may be analytically continued across 
this segment into the supersonic region of the flow. Hence, the 
flow in the subsonic region determines the flow in the supersonio 
region near the transition line. Equations (1.9) retain their 
meaning also for supersonio velocities. For, if ¥ > a, s will 
be purely imaginary and K negative. .Setting s * is and 
K <= -K in equations (1.9) gives 


&P 0 _ 

ds 


( 1 . 11 ) 


Thus, if in the solution $0,s) and ^0,8) of equations (1.9) 
determining the flow in the subsonlo region, s is set equal to is* 
and the real parts of the expressions thus obtained are taken, the 
solution of equations (1.11) determining the flow of the gas in the 
supersonio region near the transition line is obtained. 

2. Investigation of equations (1.9) near the transonic line 
W ■ From equation (1.8), it follows r tlhat in' the plane of the 

variables 0,s the upper half -plane s > 0 will correspond to 
the region of subsonic velooity and the axis of the abscissa 
s b 0 to the line of sound velocity. It is therefore necessary 
to consider the behavior of the funotion K(s) for small values of 
the variable s. 

Equation (1.2) is represented in the form 

y2 = (H+l)a» 8 _M 2 _ 

2+(K-l) M 2 


( 2 . 1 ) 
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Substituting this value of W in equation (1.8) giveB 

B .. H — tei-i )?. 2 , it (t - aCF, h 2 .sjiN 

Jo a-t 2 )(h 2 -t 2 ) V v *-y 


(l-t a )(h 6 -t‘) 
Computing the integral gives 


1 . /h-tX 11 1+t 

" 2 1 s [yi+v 1-t 


r 

i 


If the equation is expanded in a power series in, t. 


B.^i' t 3 + sii t 5 + £!4 t 7 + . . . 


3h 


51r 


7h 


6 


Then, 


( 2 . 2 ) 


t * a-jS J 


^ + ags 3 / 3 + a,s 5 / 3 + . . . 


a 5 s 


(2.3) 


where 



Further, from equation (1.2) and the adiabatic condition after 
simple transformations, 



Substituting this value in expression (1.10) gives for K(s) 



(2.5) 
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Substituting in equation (2.5) the series for t (equation 2.3) 
givss 

\jK(a) m b^ 1 / 3 + b 3 s 3 / 3 + b 5 s 5 / 3 

( 2 . 6 ) 

It follows from equation (2.6) that for a Velocity near that of 
sound, that is, for small s, consideration may be restricted in 

series (2.6) to the first term by setting VTT » b^e 1 / 3 . The 
fundamental equations (1.9) then assume the form 



- bjS 


1/3 d\| / 


, 1/3 3 * 


whence for the functions ty(e,s) and <p(0,s) there is obtained 


+ & + jl . o 


3s 


*!| + »ff - i *? = 0 (2.e) 

3s 2 3s 2 38 Ss 


Haying determined from equations (2.8) the velocity potential 
<p(0,e) and the stream function ty(9,s), by equations (1.5) the 
coordinates x and y are found in the flow plane. It is easily 
seen that equations (1.5) can be represented in the form 



(2.9) 


If, however, we pass from the ezact equations (1.9) to the approxi- 
mate equations (2.7), the expressions (2.9) cease to be exact dif- 
ferentials. Hence, simultaneously with the passing from equa- 
tions (1.9) to equations (2.7) it is necessary to introduce 


1184 



HACA !3M 1236 


7 


instead of the relation (1.8) between s and V a new relation 
such that the expressions (2.9) remain exact differentials. In 
order to obtain this relation in equations (2.9) are substituted, the 
values of the derivatives of the function op from, equations (2.7) 
after which the condition that dx from equations (2.9) is a 
total differential assumes the form 


a 



e . !° sin 
be p 


e 


aV\ i 
ds yw 


a 

3s 



e + £2 sin a 

ds p be J w 


Carrying out the differentiation and making use of the first of 
equations (2.8) satisfies this condition if the following equations 
are satisfied (the expression for dy likewise then becomes a 
total differential) : 


fo 

pW 



** w JL fPo) 

V dB [pvj 


Setting s 


2 3/2 

tj gives 






2 \l/= V 4 fo 

w = PW 


(2.10) 


Thus, the equation determining l/W 


is obtained, namely, 



( 2 . 11 ) 


The functions satisfying this equation are called Airy functions. 
Tables of these functions have been computed by 7. A. Pock 
(reference 5) . Thus 


| = C ± u(q) + C 2 v(tj) 


( 2 . 12 ) 
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where u(tj) and v(q) are two linearly independent tabulated 
integrals of the equation (2.11). The constants of integration are 
determined from the conditions 



± = /3\ 1/3 p 0 y 

dT ) Wi)«0 W 'b 1 p # a # “ a* 



where the latter relation is obtained from equations (2.10) and 
(2.6). After computation, 


C 1 = [n(0) - v'(0)| C 2 - i [u ' (0) - yu(0)j (2.13) 


where from reference 5, 


u(O) + iv(o) = -LhK 


3 4y/3 r(2/3) 


- in 
exp — 


u ' (0) + iv'(O) = — Lf* exp —■ 

T 75 ^ 6 


From the first of equations (2.10)' and (2.12), 

p 0 AN 1 / 3 C u'(tj) + Cj'(n) 

■Z— *= ( 3-) b l — T\ (2.14) 

P \ 3 / C^uOl) + C 2 v(q) 

The expressions (2.12) and (2.14) thus found for 1/fy and p Q /p 
must be substituted in equations (2.9). 

3. Determination of flow in feed part of de Laval nozzle. - 
The shape of the de Laval nozzle was determined such that the 
distribution of the velocity over the section tended, with increasing 
distance from the critical section upstream of the flow, to a uni- 
form flow with a certain subsonic velocity Wq near the velocity 



NA.CA 054 1236 


9 


of sound. It was necessary that the walls of the nozzle for 
x-> - » have a horizontal asymptote y = ±H (fig. l) , the mag- 
nitude H "being determined from the amount of gas flow and the 
velocity Wq . 


In the plane of the variables 0,s, there corresponds to the 
infinitely distant section in whioh the velocity is constant and 
equal to Wq, the point A with coordinates 6=0, s = Sq. To 

the lines of flow there correspond in' the 6,s plane a bundle of 
curves issuing from the point A. In order to obtain a solution 
of equations (2.8) having the stated properties, in the upper half- 
plane of 6,s bipolar coordinates are introduced (fig. 2). 


a 


log 


e 2 +(s+s Q ) 

e 2 +(s-s 0 ) 2 


P = arc tg — — — — £ (3.1) 

e d +s -Sq 


Thus the lines a = constant constitute a family of circles 
with centers on the s-axis and the lines (3 = constant constitute 
a family of circles with centers on the 6-axis passing through the 
point A (fig. 2). The equations of the families of these circles 
are, respectively, 


g G 2 

0 2 +(s+s n cth a) 2 = — (6-s n ctg p) 2 + s 2 = — — 
u sra u sin 2 P 

The first of equations (2.8) transformed into bipolar coordinates 
has the form 


a 2 ^ + a^ + i 

2 2 3 sh a(ch a+cos P) 

oa op 


, . d'if 

(1+ch a cos Pjv- 1 + sh a sin P - 


= (ch a + cos p) 
Equation (3.2) is reduced to the form 


1/6 


da 


X(a,p) 


ft 

jj 


dp, 

(3.2) 

(3.3) 


x. 2 * 2„ + 3 ot + 

oa d p 


_1 

36 


x= o 


(3.4) 


in which the variables are separable. 
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In seeking a solution of equation (3.4) of the form 
X= X(a) Y(p), the ordinary equations for X(a) and Y(p) 
are obtained: 



+ n 2 Y = 0 


d 2 X j. cth a dX , f 1 „2\ v _ 0 

^ -x- as + {m - ) x - 0 


where n is an arbitrary constant. If n is set equal to 0 , 


d 2 Y 0 dfx cth a dX _1_ x _ Q (3.5) 

2 2 3 da 36 “ 

dp da 


The first of equations (3.5) gives Y = C, + C ? P and the second, 

O X C» 

by the substitution t = ch a, reduces to the hypergeametric 
equation 


t(l-t) 




(3.6) 


the general integral of which can be represented in the form 


X = C, 


.- 1/12 


F {b b -I) - w (A- b h *) 


If - 0 and considering equation (3.3), the required solu- 

tion of equation (3.2) is in the form 



a +• cos 
ch a 



F 





(3.7) 


Returning to the initial variables 9 and s according to 
equations (3.1), after simple transformations 



1184 


NACA TM 1236 


11 


^ n (9,s)= 


2s r 


.2 2 2 
\0 +e +8 q 


y /6 r (i Al [8^( S - eo j [a^(.->. 0 ) g ] \ 

2/ \1 2 l 2 ( 0 2 + S 2 + s o 2 ) 2 / 


tan-1- 


2s 0 e 


e 2 +a 2 -s 0 2 


(3.8) 


This solution corresponds to the flow of gas in a nozzle having the 
shape shown in figure 1. Such flow cannot, however, be continued 
into the region of supersonic velocities. In order that a certain 
subsonic flow having a straight streamline may be continued into 
the supersonic region, it is necessary and sufficient, as has been 
shown by F. I. Frankl (reference 2) (see also S. A. Christianovich, 
reference 3, ch. Y.), that the stream function ty(0,s) have on 
the transition line the form 


^(0,0) = a^ 1 ^ 3 + A 3 e 3//r3 + a 5 s 5 ^ 3 + . : . (3.9) 

The solution (3.8) does not, however, satisfy this condition inasmuch 
as on the transition line (s = 0) it has the form 


4r 0 (©,0) 



1>x ) 


tan -1 


2sq6 


0 2 



In order to continue the flow with the stream function of equa- 
tion (3.8) Into the supersonic region, It is necessary to add to 
equation (3.8) the solution of the first of equations (2.8) satis- 

l/3 

fying the condition ^(0,0) » 6 ' . Such a solution, as is shown 
in reference 4, has the form 


A 3 ^ 3 + SAt)^ - 30 = 0 


whence setting 



we obtain 



2/3 



(3.10) 
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Hence, In order to construct the flow in the de Laval nozzle 
having the shape shown in figure 1, it is sufficient for the stream 
function \}/(e,s) to assume the form 

We,a) = A 0 Aj/ 0 + A^ + A 3 ^ 3 0k 3 = 8) (3.11) 

where Aq, Aj_, and Aj are arbitrary constants and the functions 
and ^ are given by equations (3.8) and (3.10). 

Equation (3.11) for ^(0,s) is obtained if in the expansion 
(3.9) the first two terms are retained. With the values of the 
constants Aq, A^, and A^, a nozzle may be constructed sufficiently 

near the given nozzle of the shape under consideration. 

The equations determining the functions ^ and ^ in the 

supersonic region will be determined. In the expression (3.8) for 
s = 0, the argument of the hypergeometrie function attains the 
value unity and for the supersonic velocities, that is, for Imagi- 
nary s, although remaining real, becomes greater than unity. 

The formula giving the analytic continuation of the hyper- 
geometric series (reference 6) is used 

F (a,b,c,t) * F (a,b,a+b-c+l, 1-t) + 

r(c-a) r(c-b) 

(l-t) G " a- k F (c-b, c-a, c-a-b+1, 1-t) 

r(a) r(b) 

which in the case considered has the form — 

F (p 1 > *) = TTlT/T2^57l2T F 5P I’ 1_ *) + 

T (l/W/(7/l2)‘ (l " t)l/3F (p' S' t' 1-t ) 

and the characteristic coordinates X= 6 - is and p. =_0 + is 
are determined. After computations, the following equations are 
obtained : 


fQXl 
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2a, 


\l/6 




r(l/3) w /■ 1 7 2 ■s o 2 (ji->0 / 

' — ’ ■ ■» 


r(liyi2) r(5/12) \12* 12*3* " (^+3^2^ 


/Sq(^-X) 

«1/12) t*(7/12) L +s 2 




\2/3 


-ip / 5 11 4 
^ lI2 , 12 , 3 , 


3q(^-X) 2N 


2\2 


(X^+Sq) 


arc tg 




28q(X-h0 

Xn-s 0 2 

(3.12) 

(3.13) 


These expressions determine the flow in the regions 1 and 2 between 
the transition line and the characteristic passing through the center 
of the nozzle and directed upstream, of the flow (fig. 1) . The fur- 
ther computation of the supersonic part of the nozzle can be carried 
out by the method of characteristics. 


Translated by Samuel Reiss 
National Advisory Committee 
for Aeronautics. 
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